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Abstract: The lunar Cherenkov technique, which aims to detect the coherent Cherenkov radiation pro- 
duced when UHE particles interact in the lunar regolith, was first attempted with the Parkes radio- 
telescope in 1995, though the theory was not sufficiently developed at this time to calculate a limit on 
the UHE neutrino flux from the non-observation. Since then, the technique has evolved to include ex- 
periments utilising lower frequencies, wider bandwidths, and entire arrays of antenna. We develop a 
simulation to analyse the full range of experiments, and calculate the UHE neutrino flux limit from the 
Parkes experiment, including the directional dependence. Our results suggest a methodology for planning 
future observations, and demonstrate how to utilise all available information on the nature of radio pulses 
from the Moon for the detection of UHE particles. 



Introduction 

The lunar Cherenkov technique aims to detect 
ultra-high energy (UHE: >EeV) cosmic rays and 
neutrinos via the coherent Cherenkov radiation 
emitted in the radio regime by their interactions 
in the outer layers of the Moon. The radiation 
is detected via ground-based radio-telescopes, as 
first proposed by Dagkesamanskii & Zheleznykh 
[1] and attempted at Parkes [2]. 

The use of existing radio-telescopes has meant that 
the choice of observational parameters such as fre- 
quency, bandwidth, and length of run has been 
limited. A 'next-generation' radio-array, the SKA 
(Square Kilometre Array, [3]), is currently in the 
planning stages, providing a unique opportunity 
to optimise a powerful instrument for this techn- 
qiue. The ultimate goal of the lunar Cherenkov 
technique is to determine the energy and arrival 
directions of the highest energy cosmic rays and 
neutrinos. Here, we propose how observations 
might be tailored to achieve this, by examining the 
frequency- and directional-dependence of the sen- 
sitivity of both the Parkes experiment and a nomi- 
nal SKA to UHE neutrinos. 



Simulations 

The program used to simulate the lunar Cherenkov 
technique was a Monte Carlo code described in 
James et al. [4]. This code, like previous simu- 
lations, only included a single layer in which de- 
tectable UHE particle interactions were allowed. 
However, as the regolith properties are expected to 
reflect that of the underlying material, we modified 
the code to allow UHE neutrino interactions in the 
sub-regolith layer, modelled with density 3 g/cm 3 . 
This layer was treated simply as denser regolith, 
with refractive index n and absorption length I 
scaled by density as per Olhoeft & Strangway [5] 
to n = 2.5 and I = 29 A. 

The effect of including only the classical regolith 
was to artificially limit the aperture to neutrinos of 
energies sufficiently large that interactions at the 
base of this (typically 10 m deep) layer were be- 
ing readily detected, especially at low frequencies, 
where attenuation within the regolith is minimal. 
Including two layers is more appropriate than ex- 
tending a single layer to greater depths, since some 
radiation losses are to be expected during the trans- 
mission of radio-waves from the underlying mate- 
rial (be it lava flows in the mare or megaregolith in 
the highlands) to the regolith. 
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The parameters used to simulate the Parkes exper- 
iment are also described James et al. For the nom- 
inal SKA, we use an effective collecting area of 
A c g — 1 km 2 , a system temperature T sys = 50° K, 
and uniform coverage of the Moon over all fre- 
quencies. Detection required an integrated signal 
strength in V/m over the full bandwidth of 8 times 
the RMS noise voltage. 

Results - Frequency Dependence 

The first lunar Cherenkov experiments observed at 
frequencies in the > 1 GHz range, as the power 
radiated at the Cherenkov angle scales with v 2 un- 
til the turnover frequency (2.3 GHz in the regolith 
[6]), above which decoherence effects from within 
the particle cascade become important. Recently, 
focus has been on observations at lower frequen- 
cies, as proposed by Scholten et al. [7]. At low fre- 
quencies, the Cherenkov cone is broader, and the 
radiation suffers less attenuation in propagating to 
the surface, allowing signals from a wider range of 
interaction geometries at greater depths in the re- 
golith to be detected. For a simulated experiment 
with LOFAR, assuming a bandwidth of 20 MHz, 
Scholten et al. found an optimum frequency range 
in the 115 - 240 MHz band, below which (in the 
30 — 80 MHz band) the effect of increasing sky 
temperature on sensitivity became significant. 

High-frequency experiments have a lower energy 
detection threshold, since the signal from geomet- 
rically favourable events peaks in the GHz regime. 
How this impacts upon event rates depends on the 
shape of the UHE particle spectrum, and the aper- 
ture at these energies. Another reason to not limit 
observations to low frequencies only is that to per- 
form true UHE particle astronomy, the type, en- 
ergy, and arrival direction of any detected particles 
will need to be determined. The sensitivity of the 
high frequency component to interaction geome- 
try means it could be used to resolve ambiguities 
between interaction energy, interaction depth, and 
orientation of the shower axis with respect to the 
observer, which would likely remain unresolved in 
an experiment with a low frequency range. 

The range of spectra plotted in Fig. 1 shows the 
dilemma. Any single choice of bandwidth over 
which to trigger results in reduced sensitivity to 
some fraction of events, since noise from the fre- 
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Figure 1: Simulated spectra E(y) of four different par- 
ticle interactions as viewed from Earth. The three pa- 
rameters (E, d, 9) give interaction energy, depth (m), 
and viewing angle from the Cherenkov angle respec- 
tively. 

quency range over which the signal is weakest will 
force a higher detection threshold. To detect the 
full range of events, it will be essential to observe 
across the entire frequency range. However, the 
design of an appropriate trigger algorithm is non- 
trivial, and as discussed by McFadden et al. [8], is 
a work in progress. Here we restrict our results to 
'small' (At/ = 0.3 v) bandwidths, which necessar- 
ily underestimates the SKA's potential, but gives a 
good indication of the relative importance of the 
frequency ranges examined. 

We run our simulation for the range of frequencies 
from 100 MHz to 2 GHz, using bandwidths set at 
a fixed fraction (30%) of the central frequency, re- 
flecting that instruments designed to observed at 
higher frequencies tend to have a higher available 
bandwidth. The results, together with the recalcu- 
lated Parkes aperture, are shown in Fig. 2. 

For neutrino energies > 10 21 eV (i.e. well above 
threshold), and frequencies between 140 MHz and 
1 GHz, the total effective aperture varies approxi- 
mately as , in agreement with that of Scholten 
et al. However, this dependence weakens at lower 
energies, and as neutrino energy decreases below 
10 x the effective detection threshold, the aperture 
peaks at higher and higher frequencies. Therefore, 
the choice of an optimal aperture depends on the 
assumed UHE neutrino flux, about which little is 
known. 
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Figure 2: Calculated aperture to UHE neutrinos of a 
nominal SKA observing at various frequencies. Also 
shown is the Parkes aperture both with (this paper) and 
without ([4]) sub-regolith layers included. 



Results - Directional Sensitivity 

Including geometrical tags in the simulation allows 
the sensitivity to be determined as a function of 
particle arrival direction, relative to the lunar cen- 
ter and beam pointing position. Fig. 3 shows the 
instantaneous aperture of the Parkes experiment in 
limb-pointing configuration. The total fraction of 
the sky seen at any one instant is small; the Moon's 
opacity to neutrinos at UHE creates a 'shadow' 
much larger than its apparent size (0.5°) on the 
(left) side opposite the beam pointing position. 
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Figure 3: Directional sensitivity of the Parkes experi- 
ment to 10 22 eV neutrinos in limb-pointing mode. The 
Moon is at (0, 0), beam centre at (0.5° , 0). 

We define the 'directionality', d, to be the ratio 
of peak over mean directional sensitivity, which 



measures the ratio of potential point-source sen- 
sitivity over the sensitivity to an isotropic flux. 
For Parkes in limb-pointing configuration, this is 
28; should a point-like neutrino-source be sus- 
pected/discovered, a careful choice of observation 
times and beam pointing positions would enable a 
much stronger limit to be placed. 
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Figure 4: Normalised instantaneous aperture (relative 
units) of a nominal SKA (see text) observing from 100 — 
150 MHz to 10 21 eV neutrinos. The plot is over the 
entire sky, with the Moon at the centre. 

For experiments (eg. with the SKA) which observe 
the Moon uniformly, the sensitivity forms an an- 
nulus around the Moon, as shown in Fig. 4. Simi- 
lar results would be obtained by placing multiple 
smaller beams around the limb of the Moon, as 
might be achieved with a focal plane array. Unlike 
the total aperture, the directional properties vary 
slowly with observation frequency, with the annu- 
lar sensitivity increasing in width with reduced fre- 
quency. An important result is that for energies 
well above threshold, the aperture to all arrival di- 
rections is greater at low frequencies. 

Fig. 5 plots the long-term directional sensitivity to 
10 21 eV neutrinos of a low- to mid-frequency SKA 
(350 - 490 MHz) placed at the equator. Evidently, 
there exists a significant overlap between the nom- 
inal ANITA and SKA apertures, and large regions 
of the sky to which neither experiment is sensitive. 
SKA coverage will be more uniform at lower fre- 
quencies. Fortuitously, the sinusoidal 'wobble' of 
the lunar orbit about (5 = 0° has phase aligned 
approximately with the supergalactic plane, from 
which we might expect to see an excess of UHE 
particles. Since the SKA will be sited in either 
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Figure 5: Relative directional sensitivity of a 350 — 
490 MHz SKA located at the equator to 1 ZeV neutri- 
nos (contour increment 0.1). Observations were spread 
uniformly over a lunar month. Also shown is the peak 
ANITA sensitivity range (~ 99% within -10° < S < 
+15° [9]), the galactic (solid white line) and super- 
galactic (dotten black line) planes, the apparent lunar 
orbit (thick black line), and the positions of some no- 
table AGN, used as proxies for point-like UHE neutrino 
sources. 

South Africa or Australia, the sensitivity will be 
greater to any Southern Hemisphere sources (eg 
the Galactic Centre), since the Moon will be vis- 
ible for longer when it is in negative declinations. 

Discussion 

The highly directional sensitivity of the Parkes ex- 
periments suggests that current limits on the UHE 
neutrino flux are a strong function of celestial 
coordinates. The current strongest limits on an 
isotropic flux below 10 23 eV come from ANITA- 
lite [10], which had little sensitivity beyond 15° 
from zero declination. Limits from other lunar 
Cherenkov experiments observing at higher fre- 
quencies at Goldstone [11] and Kalyazin [12] will 
be more directional than those at Parkes, and thus 
large areas of the sky would still be left unprobed. 
This leaves ample scope for experiments with dif- 
ferent directional sensitivities to make useful ob- 
servations to further constrain the UHE neutrino 
flux. 

The future of the lunar Cherenkov technique as a 
means to determine the flux of the highest energy 
cosmic rays and neutrinos hinges on the ability to 
turn positive detections into useful information on 
particle energy and arrival direction. This requires 
taking full advantage of the wide frequency range 



available to instruments such as the SKA, and will 
likely require multiple detection algorithms, opti- 
mised for various pulse spectra, to be used in par- 
allel. The optimisation must take into account the 
steepness and degree of anisotropy in the (known 
or expected) UHE particle flux, as well as the sen- 
sitivites of other experiments. The already high 
signal processing requirements of de-dispersion 
and sub-nanosecond time resolution over a large 
array will make further innovation in pulse detec- 
tion technology a necessity. 

Conclusion 

We have found the lunar Cherenkov technique to 
be highly directionally sensitive, especially at high 
frequencies. We have also demonstrated that ob- 
servations over a broad frequency range with the 
SKA will be a very powerful technique for UHE 
particle astronomy. With further improvements 
in simulation methods, and knowledge of the lu- 
nar near-surface, it will be possible to devise ad- 
vanced pulse detection and reconstruction algo- 
rithms, based on the methodology presented here. 
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